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ABSTRACT 


A  software  system  for  a  PDP-7  digital  computer  with  a  cathode  ray  tube  display 
has  been  designed  to  process  seismic  data.  The  system  permits  quick  visual  inspection 
of  digitized  data  and  allows  easy  application  of  powerful  programs  which  operate  on  the 
digitized  data  or  on  the  results  of  previously  used  programs.  Some  operations  which 
can  be  performed  are:  epicenter  location,beamforming,  magnitude,  complexity  and 
spectral  ratio  computation,  filtering,  autocorrelation,  Fourier  transformation,  sono¬ 
gram  generation  and  automatic  event  detection.  A  human  operator  is  in  the  processing 
loop,  inspecting  the  output  at  each  step  before  applying  the  next.  This  system  has 
greatly  increased  the  speed  and  efficiency  of  much  of  our  seismic  data  processing. 
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I.  GENERAL 


There  is  a  trend  for  data  acquisition  and  other  field  systems  to  use  digital 
magnetic  tape  as  the  recording  medium  in  lieu  of  analog  devices  such  as  chart,  film  or 
analog  tape  recorders.  This  trend  is  likely  to  increase  as  cheaper  and  more  reliable 
digital  tape  recorders  become  available  and  processing  by  digital  computers  becomes 
more  popular.  Digital  tape  recording  has  some  large  advantages  over  most  analog  re¬ 
cordings,  viz.  multichannel  capability,  precise  timing,  increased  dynamic  range, 
improved  accuracy,  *  and  compatibility  with  digital  computer  input  devices.  Experi¬ 
menters  are  reluctant  to  use  digital  systems,  however,  because  of  the  increased  cost 
and  complexity  concommitant  with  these  systems  and  also  because  it  is  difficult  to 
visually  examine  the  data.  This  paper  describes  a  system  which  not  only  circumvents 
the  latter  difficulty,  but  makes  standard  data  processing  schemes  that  require  the  full 
power  of  a  digital  computer  easy  to  use,  apply,  and  interpret. 

Toe  Lincoln  Data  Analysis  Console  is  a  computer -oriented  system  which  is 
designed  to  analyze  seismic  data  from  the  Large  Aperture  Seismic  Array  (LASA).  * 

(It  can  be  used  to  analyze  any  data  which  has  been  digitized  onto  magnetic  tape  in  any 
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one  of  the  several  LASA  formats.  )  This  system  is  designed  to  do  many  of  the  standard 
data  processing  tasks  more  efficiently  and  accurately,  and  yet  be  quicker  and  easier  to 

*  Off-the-shelf  digital  systems  claim  over  80  db  dynamic  range  with  0.01%  accuracy, 
while  it  is  difficult  to  get  more  than  60  db  and  0.  1%,  respectively,  from  analog  recordings. 
It  is  interesting  to  note  that  the  dynamic  range  and  accuracy  of  digital  systems  is  not 
limited  by  the  recording  method  per  se,  but  by  the  anulog-to-digital  converting  processes. 
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use  than  a  "batch  processing"  computer  of  the  type  that  is  found  in  most  modern  data 
centers.  The  system  uses  a  cathode  ray  tube  to  display  the  input  data  or  the  computer 
output  quickly  and  in  an  easily  interpretable  manner  to  a  human  data  am  st  who  can 
alte’*  the  data  or  give  new  commands  to  the  computer  via  a  fiber  optics  light  pen,  ad¬ 
justment  knobs,  or  a  teletypewriter  keyboard  (see  Fig.  1). 

This  man-machine  system  puts  a  man  in  the  feedback  loop  to  do  pattern 
recognition  and  to  make  decisions  based  on  his  experience  as  a  trained  seismologist 
or  geophysicist,  and  to  do  other  tasks  that  are  very  difficult  for  a  modern  electronic 
computer.  The  computer  is  in  the  loop  at  a  place  where  it  can  function  with  an 
efficiency  equal  to  or  better  than  that  of  a  human. 

One  of  the  difficulties  of  processing  geophysical  data  in  large  computing  centers 
arises  from  the  large  amounts  of  both  input  data  that  is  required  and  output  data  that  is 
produced.  With  data  from  a  large  seismic  array  such  as  the  Montana  LASA,  the  input 
might  typically  be  a  reel  of  magnetic  tape  and  the  physical  output  might  be  25  feet  of 
32 -channel  chart  recording,  of  which  only  25  constants  are  useful  output  (which  may  be 
the  input,  along  with  the  original  data  tape,  to  the  next  step  in  the  processing).  Each 
iteration  or  pass  through  the  computing  center  takes  time,  generates  a  lot  of  superfluous 
data,  and  can  introduce  error  in  each  two-way  digital -to -analog  and  analog -to -digital 
conversion  process  (e.g.  errors  in  reading  the  chart  recorder  and  transcribing  to 
punched  cards). 
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The  Lincoln  Data  Analysis  Console  is  designed  tc  eliminate  these  difficulties  by 
presenting  to  the  operator  quickly  and  simply  only  that  data  which  he  desires  to  see. 
This  data  is  presented  in  a  quasi-analog  manner,  i.e. ,  it  appears  analog  to  the  eye 
and  can  be  manipulated  by  analog  devices  (e.g.  light  pen  and  knobs),  yet  the  full  ac¬ 
curacy  of  the  digital  values  is  retained  for  computation  because  the  data  never  leaves 
the  computer.  For  example,  the  operator  can  easily  scan  the  input  data  tape  to  insure 
it  is  not  only  good,  but  it  is  indeed  the  data  he  desires  to  have  processed.  Useless  data 
is  not  put  onto  paper,  but  it  is  readily  available  for  perusal  if  need  be.  Only  after  the 
data  has  been  reduced,  or  after  the  operator  decides  he  really  wants  it,  is  a  permanent 
hard  C077  produced  in  the  form  of  chart  recordings, 'scope  photographs,  teletype  copy 
or  punched  paper  tape. 

The  system  is  organized  so  that  the  operator  can  have  all  the  processing 
programs  that  have  been  written  and  debugged  available  literally  at  his  fingertips. 

This  ability  to  process  the  results  of  a  prior  processing  without  punching  data  cards 
and  submitting  for  another  pass  through  the  computing  center  saves  time.  Typically, 
operations  that  take  several  days  via  other  methods  can  be  done  in  a  few  minutes  on 
the  analysis  console. 
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II.  THE  '"WTOR  SYSTEM 


The  operation  and  control  of  the  console  are  done  by  a  master  program  which 
is  called  the  MONITOR  system.  This  program  always  resides  in  the  computer  core 
memory  and  controls  the  calling  and  execution  of  all  the  commie  system  programs 
which  are  stored  on  magnetic  tape  01  a  drum  as  a  program  file.  The  MONITOR  pro¬ 
gram  always  knows  which  programs  have  been  used  and  hence  which  physical  param¬ 
eters  have  already  been  determined.  The  MONITOR  saves  all  these  parameters  so 
that  further  processing  will  have  access  to  them.  Any  processing  that  tries  to  use 
parameters  which  have  not  yet  been  determined  will  not  be  executed  by  the  MONITOR 
until  they  have  been  determined.  The  MONITOR  program  is  the  only  program  that 
must  be  manually  read  into  the  computer.  This  is  done  once  at  the  start.  A  self¬ 
protect  feature  makes  it  impossible  to  b_  destroyed  by  any  of  the  software. 

All  commends  to  the  MONITOR  are  input  to  the  computer  via  the  teletypewriter 
keyboard  in  a  one -character  mnemonic  code.  The  MONITOR  will  fetch  the  proper 
program  from  the  program  file,  transfer  it  into  the  core  memory  and  execute  it.  All 
the  programs  are  arranged  so  that  they  return  to  the  MONITOR  automatically.  Control 
can  be  returned  to  the  MONITOR  prematurely  by  typing  R  (for  return)  on  the  keyboard. 
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HI.  THE  DISPLAY  SYSTEM 


The  console  display  system  is  the  basic  program  from  which  most  of  the  other 
programs  are  called.  It  is  designed  to  make  digitized  seismograms  accessible  to  a 
seismologist  in  a  quasi  analog  manner  to  which  he  is  accustomed.  This  program  dis¬ 
plays  the  seismograms  (earth  motion  vs^time)  on  the  scope  with  six  knobs  controlling 
the  parameters  of  the  display. 

The  six  kncbs  control  the  horizontal  and  vertical  gain  of  the  seismograms,  Lhe 
vertical  separation  between  the  seismograms,  the  horizontal  and  vertical  nosition  of 
the  seismograms,  and  a  clipping  level.  Figure  2  shows  some  typical  examples. 

The  operator  can  use  this  program  to  scan  the  data  visually  and  pick  a  portion 
of  any  seismogram  for  further  processing  by  other  programs.  This  picking  is  done  by 
pointing  the  light  pen  at  any  waveform  which  causes  it  to  be  saved  in  a  "reference 
buffer.  "  The  contents  of  this  reference  buffer  will  be  displayed  at  the  top  of  the  screen 
with  a  vertical  cursor  appearing  at  the  exact  point  in  the  waveform  that  was  touched  by 
the  light  pen-  This  "reference"  trace  is  unaffected  by  most  of  the  knobs.  Thus  one 
can,  for  example,  move  all  the  traces,  one  by  one,  alongside  this  reference  by  adroit 
manipulation  of  the  horizontal  and  vertical  position  knobs  to  do  visual  alignment  and 
correlation.  This  process  is  used  to  pick  arrival  times  for  the  various  traces  to  be 
used  for  beamforming  and  location  (see  below).  By  proper  adjustments  of  the  gains, 
this  time  picking  can  easily  be  done  with  a  precision  of  one  sample  (50  milliseconds 
when  using  short-period  data)  as  illustrated  in  Fig.  3. 
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IV.  THS  CONSOLE  PROGRAMS 

Some  of  the  operations  that  can  be  done  with  the  analysis  console  will  be  briefly 
described.  *  The  3^ stem  is  designed  so  that  additional  programs  can  easily  bj  added 
when  they  are  written  and  checked  out. 


A.  INITIALIZE 

This  is  usually  the  first  program  that  is  called  and  executed.  It  reads  the  data 
tape  and  determines  which  of  the  several  allowable  formats  the  tape  has  been  recorded 
in.  (Each  site  may  have  its  own  tape  format, )  The  ptogram  types  out  this  information, 
along  with  the  date  and  starting  time  of  the  data  tape.  The  operator  can  select  which 
data  channels  he  wants  to  operate  upon  (there  may  be  up  to  several  hundred  channels 
recorded  on  the  data  tapes)  and  what  sampling  interval  should  be  used.  The  sampling 
can  be  changed  by  instructing  the  computer  not  to  use  every  sample  that  was  recorded, 
but  to  skip  a  fixed  number  of  data  samples  from  each  channel.  This  program  initializes 
all  the  other  programs  to  use  these  pre-st  acted  parameters. 

B.  LOCATION 

If  three  or  more  arrival  times  are  picked,  either  by  the  display  program  or  by 
manually  typing  them  in,  the  best  fitting  plane  wave  that  minimizes  the  r.  m.  s.  error 

*  A  detailed  description  of  ail  the  programs  is  availabl  Tom  the  author  on  request. 
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will  be  calculated.  From  this,  the  horizontal  phase  velocity  and  azimuth  will  be 


calculated.  In  addition,  the  relative  arrival  times  when  tins  plane  wave  would  hit  all 
21  sites,  along  with  the  residuals  from  these  times,  is  typed  out.  Provision  is  made 
for  deleting  times  from  any  sites  that  may  have  large  residuals  and  repeating  the  cal¬ 
culations.  (If  the  data  is  from  short -period  seismometers,  station  corrections  will 
have  automatically  been  inserted  in  this  calculation,  and  the  distance  from  the  center 
of  LASA,  the  latitude  and  longitude  of  the  epicenter  will  be  computed.  Also,  G.M.T. 
time  at  which  that  portion  of  the  seismogram  corresponding  to  the  reference  cursor 
must  have  originated  will  be  calculated,  assuming  a  focal  depth  of  33  km. )  Two  con¬ 
stants  giving  the  major  and  minor  axes  of  the  location  uncertainty  ellipse  are  also 
typed  out.  Provision  is  made  for  manually  deleting  the  station  corrections.  Figure  4 
shows  the  output  for  the  event  shown  in  Figs.  2  and  3. 

C.  BEAM  FORMING 

If  one  or  more  time  picks  are  made  via  the  Display  program  or  the  manual 
type-in  program,  the  delay-and -sum  beam  will  be  formed  and  placed  in  the  reference 
buffer  as  in  Fig.  5.  From  here,  it  can  be  filtered  (see  below)  or  used  as  a  reference 
to  pick  arrival  ’dmes  from  weak  events.  Channels  that  have  been  deleted  in  the  Locate 
program  (above)  will  also  be  automatically  deleted  from  the  beam. 
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D. 


CONSTANTS 


This  program  displays  the  reference  waveform  on  the  scope.  Two  points  can 
be  selected  by  the  light  pen  and  the  G.M.T.  time  corresponding  to  each  point  will  be 
displayed  to  the  nearest  0. 1  second.  In  addition,  the  time  difference  between  the  points 
and  the  amplitude  difference  in  millimicrons  corresponding  to  the  two  points  picked 
(assuming  normal  calibration  for  the  data)will  be  displayed.  The  magnitude  ba:;ed  on  these 
last  two  quantities  will  be  typed  out  on  request  if  the  distance  has  already  been  com¬ 
puted.  Provision  is  made  for  inputting  the  distance  if  it  has  not  been  computed.  If  the 
data  is  from  a  short-period  instrument,  the  "Q"  factors  for  a  depth  of  33  km  are  used. 
See  Figure  6. 


E.  FILTERING 


The  data  stored  in  the  computer  (21  seismograms,  each  600  samples  long)  or 
the  reference  can  be  filtered.  Currently,  one  has  the  choice  of  two  filters  which  have 
been  found  useful  in  the  past  for  short -period  data. 

1.  Notch  Filter 

This  is  a  low  distortion  filter  with  two  notches  in  the  frequency 
response  (at  0.2  and  0.3  Hz)  to  attenuate  microseismic  mise 
occurring  around  these  frequencies.  Zero  frequency  is  attenu¬ 
ated  about  10  db  and  frequencies  higher  than  about  1  Hz  pass 
through  with  little  attenuation  oi  phase  shift. 

2.  Butterworth  Filter 

This  is  a  three -pole  bandpass  filter  with  3  db  points  at  0.  6  and 
2.0  Hz.  Both  zero  frequency  and  high  frequencies  (-30  db  at  5.0 
Hz)  are  attenuated,  but  some  phase  distortion  occurs. 

Figure  7  shows  some  examples  of  the  use  of  these  filters. 
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F. 


PLOT 


This  program  plots  ail  the  waveform  data  stored  in  the  computer  on  the  Sanborn 
chart  recorders  so  that  hard  copy  can  be  obtained.  This  gives  a  much  larger  scaling 
than  photographs  of  the 'scope  can. 

G.  TRANSFORM 

This  program  computes  and  displays  the  autocorrelation  function  of  the 
reference  waveform  (512  data  points),  which  has  been  modified  as  follows: 

0  j<  A,  j  >  B 

fREF  A  ^  j  s  B  j  =  1,  512 

where  A  and  8  are  set  to  any  values  via  the  light  pen. 

The  finite  Fourier  transform  of  either  the  modified  reference  or  of  its 
autocorrelation  function  (both  considered  periodic  with  a  period  of  512  samples)  can 
be  computed.  The  autocorrelation  function  can  be  modified  by  an  adjustable  rectangular 
window.  The  periodogram  of  this  transform  is  normalized  and  displayed  (db  vs  frequency) 
and  it  can  be  output  onto  paper  tape.  Figure  8  shows  some  examples. 
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H.  COMPLEXITY 

This  program  computes  and  types  out  the  complexity  oi  the  reference  trace 
(beam  or  seismogram)  starting  at  the  cursor.  The  complexity  is  defined  as  follows: 


C  = 


t  +  35  sec 
c 


2 


j  =  t  +5  sec 
c 


Vf> 


t  +5  sec 

C  2  (t-b 

j  =  tc 


where  t 

c 


is  the  time  of  the  cursor  that  was  set  by  the  light  pen,  and 


f  = 


20  . 
J 


v- 

2 


t  -20  sec 
c 


f.  . 
J 


I.  SONOGRAM 

A  sonogram  is  a  display  of  signal  intensity  (represented  by  several  different 
levels  of  brightness)  as  a  function  of  frequency  (vertically)  and  time  (horizontally). 
Figure  9b  shows  the  sonogram  of  the  test  calibrate  signal  of  Fig.  9a. 

This  sonogram  uses  a  bank  of  50  constant  bandwidth  filters,  each  with  a  (sin  x/x) 
type  frequency  response,  with  impulse  response  exactly  10  seconds  in  duration.  The 
filters  are  spaced  0. 1  Hz  apa  t,  and  they  cover  a  range  of  0.  1  to  5.0  Hz.  The  input  to 
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these  filters  is  the  reference  trace  starting  1200  samples  before  the  data  displayed  in 
the  reference  buffer  and  continuing  for  3200  samples.  The  sonogram  time  scale, 
therefore,  covers  a  range  of  160  seconds  if  short -period  data  is  used  or  over  2  hours 
if  long -period  data  is  used. 

Each  output  of  the  50  filters  is  full -wave -linear  rectified  to  form  the  signal 
intensity,  *  and  passed  through  two  low-pass  filters  in  cascade.  The  first  is  an  enerr  y 
summing  "integrate  and  dump"  filter.  It  sums  the  output  of  the  rectifier  for  20  sam¬ 
ples,  then  resets  itself  to  zero  and  repeats.  The  second  filter  takes  these  outputs, 
sampled  just  before  the  dump  and  performs  "RC"  type  low -pass  filtering.  The  3  db 
cutoff  of  this  latter  filter  is  0.  11  Hz.  The  computer  stores  the  output  of  this  bank  of 
50  "RC"  filters  in  a  two-dimensional  array  of  power  vs  frequency  and  time.  The  sono¬ 
gram  is  simply  the  pseudo  three  dimensional  display  of  this  two  dimensional  array. 

If,  for  every  time  element  in  the  sonogram  array,  the  intensity  at  every  freqv  ?ncy 
is  summed,  a  one -dimensional  "ene-gv  profile"  array  is  formed.  This  program  cal¬ 
culates  and  displays  this  energy  profile  as  a  series  of  160  fine  dots  at  the  top  of  the 
display.. 

If,  for  every  frequency  element  in  the  sonogram  array,  the  intensity  Detween 
two  set  times  is  summed,  a  one -dimensional  "frequency  spectrum"  is  formed.  This 
program  allows  the  two  times  to  be  manually  set  via  the  light  pen.  The  spectrum  is 

*  Intensity  rather  than  power  is  computed  because  it  makes  less  stringest  requirements 
on  the  dynamic  range  of  the  internal  computer  calculations. 
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displayed  as  a  series  of  50  coarsely  spaced  dots  superimposed  upon  the  energy  profile. 
Both  these  latter  displays  are  automatically  normalized  to  use  all  the  room  on  the 
scope  available  to  them. 

3  4 

The  spectral  ratio  ’  can  easily  be  calculated  by  summing  the  proper  elements 
in  the  sonogram  array  and  taking  their  ratio.  This  is  done  and  the  result  is  output  on 
the  teletypewriter. 

J.  AUTOMATIC  EVENT  DETECTION 

This  program  scans  the  data  tape  and  types  out  the  time  of  all  teleseisms 

encountered.  The  actual  mechanism  by  which  teleseisms  are  distinguished  from 

5  6 

microseism ic  and  other  types  of  interference  has  been  discussed  elsewhere.  ’  This 
operation  takes  place  much  faster  than  real  time  so  that  one  can  start  with  this  pro¬ 
gram  and  quickly  find  all  the  events  on  the  data  tape,  then  analyze  each  separate  event 
by  using  the  programs  described  earlier. 

K.  CONVOLUTIONAL  MATCHED  FILTERING 

This  program  cross  correlates  the  reference  trace  with  a  matched  filter  that  is 

7 

especially  designed  to  enhance  the  S/N  of  long -period  dispersed  surface  waves.  The 
matched  filter  is  a  sinusoid  with  a  linear  frequency  modulation  (chirp).  The  matched 
filter  has  three  parameters  that  can  be  input  to  the  program  via  the  teletypewriter. 

They  are:  the  period  at  the  start  of  the  chirp;  the  period  at  the  end  of  the  chirp;  and  the 
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length  of  the  chirp.  The  maximum  allowable  length  of  the  chirp  matched  filter  is  800 
samples  which  corresponds  to  160  to  1920  seconds  depending  on  the  sampling  chosen 


in  the  Initialize  program.  Enough  of  the  input  reference  will  be  convolved  with  this 
chirp  to  give  512  samples  of  filtered  output,  which  will  appear  in  the  reference  buffer. 
Figure  10  shows  a  typical  example  of  this  type  of  filtering. 
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Fig.  1.  The  Lincoln  Data  Analysis  Console. 
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Fig.  2a.  All  21  seismograms. 


Fig.  2b.  The  top  four  seismograms  with  greater  gain  and  separation. 
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Fig.  3a.  Seismogram  No.  4  »sed  as  reference  and  aligned  with  No.  1 


Fig.  3b.  Same  except  aligned  witi;  No.  2  and  more  gain. 
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Fig.  4.  Output  of  locate  program. 
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Fig.  5.  Beam  formed  from  all  21  seismogTams  at  top. 
Notice  improved  S/N. 


Fig.  6.  Seme  constants  for  top  seismogram  of  Fig.  2. 
Note  that  the  vertical  lines  can  be  moved  at  will  with  the 
light  pen.  The  teletypewriter  typed  MAG=5.2. 
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-64-9819 


Fig.  8a.  Top  trace  is  input  data.  At  the  bottom,  a  certain 
portion  is  selected  with  the  light  pen.  The  middle  trace  is 
the  autocorrelation  function  of  the  bottom  trace. 


Fig.  8b.  Periodogram  of  the  bottom 


trace  of  Fig.  8a. 


A 


Fig.  8c.  Same  as  Fig.  8b  showing  the  first  three  Hz  expanded. 


Fig.  8d.  Periodogram  of  the  autocorrelation  function  of  Fig.  8a 
after  it  has  been  modified  by  a  rectangular  window  (vertical  lines 
in  center  portion  of  Fig.  8a). 
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Fig.  9a.  Calibrate  signals  on  various  channels. 

-64-»»4| 


Fig.  9b.  Sonogram  of  top  trace  of  Fig.  9a.  Note  energy  profile  (fine 
dotted  line)  and  spectrum  (coarse  dotted  line)  at  the  top.  Note  slight 
third  harmonic  distortion  inadvertently  present  in  the  calibrate  signal. 
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